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We derive the essential space-time spectrum of the Maxwell equations in linear isotropic inhomo- 
geneous media together with the corresponding essential modes. These modes represent the collapse 
of the electromagnetic field into a single point in space at a single angular frequency. The location 
and frequency of the essential mode are random variables obeying the Born statistical postulate. 
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Wavefunction collapse is one of the phenomenological 
postulates of the orthodox quantum theory. However, 
this concept is so foreign to the wave behavior in general, 
that many find it impossible to accept. The recent years 
have seen an increase in the number of publications de- 
voted to the foundations of quantum physics. The more 
constructive of the new interpretations suggest that the 
collapse is just an awkward theoretical construct, claim- 
ing that there is no such thing in reality. 

We mention two such theories that seem to have passed 
the test of time and are often regarded as candidates for a 
new, less orthodox, model (see [l| for an overview). One, 
called the decoherence theory, operates almost entirely 
within the standard postulates of quantum mechanics, in- 
cluding the statistical interpretation of the wavefunction. 
The decoherence theory argues that the most important 
wave feature of the quantum behavior is the interference 
and as such it can only be observed with coherent wave- 
functions. Therefore, the apparent breakdown of wave 
behavior in the process of measurement and, more gener- 
ally, in macroscopic objects, is nothing else but the major 
loss of coherence induced by the environment. The loss 
of coherence, however, is never perfect, so that some in- 
terference is always present, even at a macroscopic level. 
Isolated microscopic objects would then never experience 
collapse at all. On the other hand, may be there is no 
such thing as a truly isolated microscopic object. 

Another theory, known as the dynamical reduction 
model, suggests that the basic wave equation of quan- 
tum mechanics - the Schrodinger equation - should be 
modified, so that it would admit collapsing solutions. 
The proposed modifications include a nonlinear and/or 
a stochastic term, since it is believed that a unitary lin- 
ear wave evolution can never exhibit a random collapsing 
behavior. The dynamical reduction program neither pro- 
vides nor aims at a perfectly collapsing wavefunction, and 
imposes the quantum statistics a priori. Hence, both the 
decoherence theory and the dynamical reduction model 
introduce something that "looks and smells like a col- 
lapse" , but in our view is not the real thing. 

We believe that the true problem with accepting the 



orthodox point of view is largely due to the absence of ad- 
equate mathematics that would describe collapse of the 
wave behavior at least in some particular setting. We 
shall not go into the details of the two interpretations 
mentioned above, which question the very existence of 
the collapse. Instead, we shall provide an example of a 
perfect collapse of wave motion, which has nothing to do 
with coherence or nonlinearity. Namely, we shall show 
that the linear Maxwell's equations admit generalized 
collapsing solutions, that are not only random, but also 
satisfy the Born statistical postulate in a certain sense. 

A few words have to be said about the mathematical 
techniques employed here. We follow up on the discovery 
of the essential spectrum in the volume integral formu- 
lation of the electromagnetic scattering problem Q de- 
rived using the symbol theory of Mikhlin and Prossdorf 
0. Our preliminary qualitative analysis indicated that 
the "essential resonance" , which happens if the essential 
spectrum gets close to zero, must be related to the so- 
called plasmon, where the electromagnetic field is coupled 
to the plasma wave and is highly localized at a plasma- 
dielectric interface This prompted further investiga- 
tion and recently we were able to prove that the mode 
corresponding to the point of essential spectrum,called 
singular mode, is a very special function, indeed jfj. We 
have used the Weyl definition of spectrum Q showing 
that the singular mode is a distribution, which is best 
described as the square root of the Dirac delta function. 
Such functions obey the Colombeau algebra || and, de- 
spite their unique properties, are not in wide use today. 

Here we apply the Weyl definition of spectrum to the 
differential Maxwell's equation and show that the elec- 
tromagnetic essential spectrum is not a unique feature of 
the volume integral formulation only. In fact, the spatial 
singular modes, derived for the integral operator in jfj, 
will appear in the differential case as well. To arrive at 
a physically interpretable picture we do not consider the 
spatial and temporal spectra separately, but immediately 
derive the essential space-time spectrum of the Maxwell's 
operator. Based on the unique structure of the essential 
modes, their natural randomness, and Born-type statis- 
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tical properties, we conclude that we deal here with the 
perfect collapse of the electromagnetic wave motion. 
Consider the Maxwell equations 



-V x H + <9 t D = J, 
V x E + d t B = 0, 



(1) 



with linear isotropic constitutive relations D(x, i) = 
e(x, t) * E(x, t), B(x, t) = /x(x, t) * H(x, f), where e(x, t) 
and /i(x, t) are some continuous functions with finite spa- 
tial support, and star denotes the temporal convolution. 
In the matrix form these equations can be written as 



(2) 
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In operator notation we simply write 
MX = Y, 



(3) 



where M is the Maxwell operator, and X and Y are six- 
vectors. We introduce the Hilbert space with the follow- 
ing norm: 



ixr dxdt. 



(4) 



According to the Weyl definition of spectrum of M one 
needs to find the sequence of functions {F n } with the 
following properties: 



lim | 

n — >oo 



IIF, 
AF r , 



(5) 
(6) 



If the Weyl sequence that satisfies these equations for 
some A has no convergent subsequence, then such A is 
in the essential spectrum. If for some other A the Weyl 
sequence does converge to a function from the Hilbert 
space, then this function is an eigenfunction, and A is an 
eigenvalue of M. 

The singular modes of the volume integral operator 
found in p] are the following vector-valued functions: 



*(or,x,x, 
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-3/4 a 5/4 



(x - x c ) exp 



(7) 



where x, x c £ K 3 and a > is the sequence parame- 
ter. The properties of these functions are summarized in 
Theorem 2.1 of The most important are the following 
two: 

f |*(a,x,x c )| 2 dx=l, (8) 

JxGR 3 

lim f /(x)|*(a,x,x c )| 2 dx = /(x c ). (9) 



This is why we refer to singular modes as the square root 
of the delta function. Further we note that the vector- 
valued \& has the form of a position vector (x — x c ) multi- 
plied by a scalar, and its Fourier transform has the form 
of an angular vector k multiplied by a scalar This 
means that the action of a curl operator on 'J gives zero 
for any value of parameters a and x c , which is the main 
reason why singular modes are possible with the first- 
order Maxwell's equation, but not with the second-order 
wave equations in homogeneous media, where the spatial 
operator is the Laplacian. 

Here we intend to take into account the temporal es- 
sential spectrum as well. Consider a similar sequence of 
scalar functions: 



<j)(/3,t,uj c ) = 

i\/27r- 1/4 /3" 3/4 texp ~ iuJ 

The Fourier transform of these functions is 

^-1/4^3/4^ _ ^ exp - L0 C \ 2 



(10) 



(11) 



and, as can be verified using the integration techniques 
developed in 0], it is normalized and also represents the 
square root of the Dirac delta function, but now in lu- 
domain, i.e., 



duj = 1, 



lim 



g(w) <j>(/3,L),u)c) duj = g(ui c ) 



(12) 
(13) 



The product of functions Q and 1|10|) can be used to 
generate the essential modes of the complete Maxwell's 
operator, i.e., 











(14) 



Substituting F c and F m in JBJ) we arrive at 



lim ||MF e - A°F°|r = |A c + ^ c £(x c ,w c )r, (15) 

Ct,0 — >-OG 



lim j 

a, (3— >oo 



A m F m|j2 = | A m + ^ c/t(XcjWc) |2 ; (lg) 



meaning that — icj c e(x c , u) c ) and — icj c /i(x c , w c ), where 
u! c E M and x c £ R 3 , are in the essential spectrum of 
the Maxwell operator. 

Obviously, the Weyl definition of spectrum allows 
working with a very broad class of functions, and it may 
seem that the modes of essential spectrum do not have 
any mathematical and as a result any physical meaning. 
Indeed, we show in 0, that the singular modes do not 
(strongly) converge to any function of the Hilbert space. 
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On the other hand, the square of these functions gener- 
ates the generalized Dirac's delta function, which in its 
turn acquires some meaning only upon integration. All 
this amounts to saying that the essential modes may al- 
low for some physical interpretation only if we square and 
integrate them. 

The physical quantity obtained by squaring and inte- 
grating the electromagnetic field over space and time has 
the dimensions of action [Js], the electric-field part of 
which is: 



S e = 



[ [ E(x,i) • D(x,i)dxdt 



— oo ./x 



Re[e(x,w)] E(x,w) 



(17) 



dxdw. 



Here we have used the power (Plancherel) theorem of 
Fourier analysis, the fact that E(x, t) is a real function 
of time, and the causality-driven assumption that the 
real part of e(x, ui) is an even function of u>, while the 
imaginary part is odd. Of course, S c is a slightly artifi- 
cial quantity. It only serves to illustrate our point and 
techniques. 

It is clear that the same numerical value of action may 
correspond to different states of the electric field. For 
example, field E(x, u) may be a "proper" solution of the 
Maxwell equations. Alternatively, it can be an essen- 
tial mode, which satisfies Maxwell's equation only in the 
Weyl sense, but produces a perfectly well-defined action. 
Consider a hypothetical transition between such states. 

Let the action Sf in the first state be represented by an 
integral over the electric field Ei satisfying the Maxwell 
equation in the usual sense, i.e. Sf is given by (|17f) with 
Ei instead of E. In the second state the action is repre- 
sented by a single essential mode 



= Co lim 

a, [3 — >oo J _ 

(j>(P,u), w c )*(a,x,x c ) 
= C Re[e(x c ,cj c )] . 



Re [e(x, lo)] x 



dxdw 



(18) 



Although, we can always make 5° = S| by picking a 
suitable value for the dimensional coefficient Co, there 
is an important difference between the "proper" field Ei 
and the essential modes. When we say that Ei satisfies 
Maxwell's equations we mean that it is a solution of a 
boundary value problem, which exists and is unique in 
most circumstances. The essential modes, however, are 
not unique. For example, for a fixed Co the conservation 
of action holds with any x c and u c at the equipotential 
hypersurface (hypervolume) of Re [e(x c , u c )]. 

Hence, the action-preserving transformation between 
the usual solutions of the Maxwell equations and the es- 
sential modes is nonunique in two respects. First, with a 
fixed value of action there are no means to decide whether 



this transformation takes place at all. Second, if such 
transformation has happened, then the parameters of the 
location of the spatial singular mode x c and the angular 
frequency uj c of the temporal essential mode are com- 
pletely arbitrary. 

Let us see, if the physical interpretation of the essential 
modes can be extended any further than the squared and 
integrated action-like expressions. The previous consid- 
eration suggests that in the transition between the states 
both x c and u c should be treated as random variables. 
Suppose that the transition is such that each essential 
mode can carry only a finite portion of action, say a, and, 
therefore, the second state corresponding to some suffi- 
ciently large numerical value of action S% ^> a consists 
of a large number of essential modes rather than a single 
mode. One possibility is to assign Co = h/so in ltTH|l . 
describing the single-mode mode transformation, where 
h is the Plank constant. Then, the amount of action dS% 
corresponding to the essential modes with positions x c 
in the subvolume dx c and frequencies oj c in the interval 
du> r is 



dS| = C/(x c )dx c #(w c )dw c x 



lim Re [e(x, u>)} x 

2 

0(/3,w,w c )*(a,x,x c ) dxdcj, 



(19) 



where /(x c )dx c and g(uj c )duj c are the relative numbers of 
essential modes with parameters within the correspond- 
ing subvolume and interval. Note that x c and uj c are still 
random. Hence, /(x c ) and g(oJ c ) can be interpreted as 
probability densities. The complete action in the second 
state will look as follows: 



S%=C lim 



f(x c )g(cj c )x 



Re [e(x, to)] x 



(20) 



0(/3,w,a; c )*(a,x,x c ) 



dx dtu dx r duj r 



Now, equating Sf — 5| again, interchanging the order 
of integration, and using the symmetry and the sifting 
property of the essential modes we obtain 



— oo J X 



c 



— oo J X 



lim 

q,/3^oo J_ 00 7 xg 



Re[e(x, to)] Ei(x, lj) 
Re [e(x, lj)] x 
f(x c )g(u c )x 



dx duj 



(21) 



(^(/3,w,o; c )*(q!,x,x c ) 



dx r duj r dx duj 



C 



Re [e(x, cu)] f(x)g(ui) dx dw. 
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We also conclude that this equality holds if 



/(x) 5 (o;) = (T 1 Ejfow 



(22) 



In other words, the probability density function of the 
random variables x c and to c is equal to the normalized 
squared amplitude of the electromagnetic field. This, 
obviously, resembles another postulate of the quantum 
theory, known as the Born rule. Although, we have to 
point out that equation (|21J) may hold for other densities 
as well, and more information is required to make 122[) 
the unique choice. 

Thus, Maxwell's equations admit generalized Weyl's 
solutions which exhibit many of the features of the 
enigmatic wavef unction collapse. Furthermore, the ex- 
istence of essential modes suggests the following hier- 
archy of equations in the electromagnetic field theory. 
There are global equations on the level of the electro- 
magnetic action allowing the hypothetical transitions 
between "proper" and "essential" states. These equa- 
tions can be satisfied both by the usual solutions of the 
Maxwell equations and by the essential modes. Mathe- 
matically, such equations are simply equalities between 
quadratic functionals defined on functional spaces far 
broader than the common and rather restrictive complete 
Hilbert space. The more informative local (pointwise) 
equations involving intensities (squared amplitudes) of 
the field are satisfied by the "proper" solutions in the 
usual sense, and only statistically by the essential modes. 
The random parameters are the locations of the spatial 
singular modes and the angular frequency of the tempo- 
ral modes. The natural emergence of the Born statistical 
rule here is truly remarkable, although, more work is re- 
quired to justify it. Finally, we have the electromagnetic 
field, which satisfies the Maxwell equations in the usual 
sense. Essential modes do not have any mathematical 
and, hence, physical meaning on that level. 

Instead of the electric field action l|17[) we could use 
a sum of several action-like or other terms quadratic in 
fields, requiring the concervation of the total expression 
under the transformation of modes. Although the exact 
physical interpretation of the transformation may change 
in that case, the methods developed here should remain 
valid. Hence, we do not know yet how to induce or control 
the electromagnetic collapse in practice, and if this can 
be done at all. We do know the condition for the essential 
resonance, i.e., A c = — iuj c e(x c , u c ) = 0. With "proper" 
eigenmodes the presence of a zero eigenvalue means that 
all energy is stored in the corresponding mode, which is 
a realizable spatial distribution of the field. The essential 



resonance, however, is not so easy to interpret. We can 
only discuss the relevant physical conditions, keeping in 
mind that x c and lo c are random. One straightforward 
case is lu c — 0, while e(x c ,w c ) ^ 0, corresponding to the 
static breakdown of wave motion. Another possibility is 
e(x c ,cij c ) = 0, while lu c =/= 0. This does not happen in 
causal media. On the other hand, with strong anoma- 
lous dispersion the real part of e may be zero and even 
negative. However, this is always accompanied by the in- 
crease in the imaginary part of e. Which rises a question: 
should not we see the absorption of the electromagnetic 
radiation, empirically described by the imaginary part 
of e, as resulting from the excitation of essential modes? 
The concervation of energy (Poynting theorem) may help 
to derive the corresponding relation. 

Hence, collapse of the wave behavior is not such a 
crazy thing after all. Mathematically it can be viewed as 
the action-preserving transition between the states cor- 
responding to discrete and essential spectra of operators, 
which does not depend on the loss of coherence due to en- 
vironment or on the presense of nonlinearity. Certainly, it 
is not clear yet how to extend these ideas on the quantum 
wave motion. In any case, our experience with electro- 
magnetics shows that the critical analysis of the mathe- 
matical conditions on the existence of the solution with 
the subsequent explicit derivation of the essential modes 
as the Weyl solutions of the first-order Dirac's equation 
may be a good starting point. 
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